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ABSTiiACT 



The positive colusnn of an Argon glow discharge was observed through a 
Fabry-Perot interferometer with a prism spectrograph in order to determine 
if a Doppler shift occurs. No Doppler shifts were observed; however, fur- 
ther refinements in equipment and technique are required before it can be 
reported that no Doppler shifts exists. 

Feasibility of removable end plate construction of discharge tubes was 
investigated and future design parameters were specified. 

Methods for obtaining emitter radial distributions from spectroscopic 
intensity measurements were investigated. In lieu of experimental data, 
artifical intensity measurements were used as inputs to the Control Data 
Corporation 1604 digital computer in order to compare and to determine the 
sensitivity of these methods. 
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1 . Introduction 
1 . 1 History 

During the past century investigations have been made of the field 
we currently call plasma physics. It was first termed the physics of gas 
discharges, then with the introduction of electronics it was called gas- 
eous electronics. The term plasma was coined in 1920 by Irving Langmuir^* ^ 
to denote a gas discharge in which quasi-neutrality is maintained, and in 
which a visible glow is apparent. 

Unfortunately, there is still no standard, all inclusive, definition 
for the term plasma. It may be defined as a neutral collection of electrons 
and positive ions; however this is an over-simplification. The true def- 
inition must be considered with respect to the application of the plasma 
and the environment under which it exists. Under these restrictions, 
the plasma may be defined in terns of Debye length, a characteristic 
linear dimension; electron and/or density temperature, or degree of ion- 
ization. 

Much of our knowledge of plasmas has been obtained using spectro- 
scopic diagnostic techniques. Using these methods, information about the 
plasma constituents, temperatures, densities, and particle velocities 
has been obtained. Recently, microwave, probe, and other non-spec tro- 
scopic methods have been perfected. Under seme circumstances they pro- 
vide more accurate data than optical devices, but they have one common 
fault, they perturb the plasma. 

As a result of the discovery of moving striation in 1874 by Wullner' 
in glow discharges, many studies of this phenomenon have been undertaken. 
Moving striations are "the alternate bands of darkness and luminosity 
traveling through the positive column of gas glow discharges in Mercury 

h/ 

and inert gases." ' ' 
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1.2 Glow Discharges. 

A large number of studies have been made of glow discharges; as 
a result, the characteristic regions of this type gas discharge have 
been classified according to the standard nomenclature shown in figure 

la. 

A glow discharge which is maintained by electrons produced at the cathode 
by a combination of positive-ion bombardment and electron emission from 
the cathode filament is called a hot cathode discharge. Current within 
the discharge is maintained by placing a potential of several hundred volts 
across the discharge tube. 

The positive column is the region of our interest. It is the ion- 
ized region that extends from the Faraday dark space to the anode glow 
region. From figures 1b, c, d, it is obvious that this region is 
characterized by high light intensity, a low potential drop, and low 
electric field strength. The visible glow from this region comes from 
transitions of bound electrons in the Argon atoms, or from free electron 
entering or leaving bound atomic states. There is no evidence that 
similar interactions are occurring with the Argon positive ions.^^ The 
physically obvious characteristics of the positive column are heavily 
dependent on many parameters such as tube geometry, gas pressure, dis- 
charge current, supplied voltage, cathode temperature, and type of gas. 

1 .3 Previous Experimental Work. 

Far too many experiments have been conducted on moving striations 

during this century to be reviewed in this paper. An excellent summary 

h/ 

of early works was made by J. J. Thompson ; the work done in the period 
J920-1956 is outlined by A. W. Cooper^^ and G. Francis^^. However, two 
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GLOW DISCHARGE CHARACTERISTICS 
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recent works will be summarized. 

Ini 

Robertson and Hakeem' ' investigated the dependence of striations 
on the presence of me castable atoms. Moving striations can originate 
either from instabilities of the plasma, or from oscillations formed at 
the plasma boundary and propagated through the plasma. Robertson and 
Hakeem concluded that in both cases, the causes for moving striations 
are realized only when metastable atoms are present in the glow discharge. 
If there are no metastable mechanisms available, such as is the case in 
the alkali glow discharges, then the instabilities are attenuated and 
moving striations are not observed. 

- /g/ 

Ascoli-Bartoli, De Angelis and Martellucci' ' used optical inter- 
ferometry methods to determine electron density distribution in a plasma. 
Their results are still meager, but the application of S ;hlieren photo- 
graphy to investigate moving striations may prove to be a very valuable 
tool in determining the phenomena that accompany them. For example, 
the analysis of density gradients, if they exist, would lead to inform- 
ation regarding the plasma flow parameters that are associated with the 
striations. 

1 .4 Theory 

Moving striations travel down thelength of a plasma column with a 

/ 3/ 

finite velocity. Cooper' ' has shown that the velocity tends to have 
several peak values over a range of pressures and currents. The stri- 
ation velocity is equal to the product of the frequency and wave length 
of the striations. 

If some of the light emitters are moving along the tube with the 
same velocity as the moving striations, a Doppler shift in the emitted 
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light should exist. In addition to this source of wave length shift, 
Doppler broadening and Stark effect cause spectrum perturbations. 

Doppler broadening is dependent on the random velocity distributions of 
the emitters within the glow discharge. This results in a smear of 
spectral lines appearing in place of a narrow, well-resolved single 
line, The Doppler shift due to the velocity of the moving striations 
should appear as a v/ell defined shift, and separation of the two 
Doppler effects is possible provided sufficiently sensitive optical 
equipments are employed. The Stark effect should not affect the result 
of this experiment, since it is related to the electric field strength, 
and in the positive column the electric field strength is low. (See 
figure Id.) 

If an emitter at rest produces a spectral line of wavelength , 
then it will produce a slightly different wave length /Y when in 
motion with velocity V. 

x\ c A ( i + it) (1.41) 

The Doppler shift ^A = A~A that results from velocity V is 
From this one sees that the velocity is 

nr ~ (1.43) 

With an anticipated emitter velocity of 90 m/s , one can expect 

o _o ° 

to observe the 8 II 5 A Argon line to shift 2.45 x. lO - "' A. In order to 
detect this magnitude of wave length shift, equipment with a resolving 
power of 3 x 10 must be employed. 
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2. Equipment and Procedures. 

2.1 Vacuum System. 

A schematic diagram of the vacuum system is shown in figure 2. 

The system was generally operated at a base pressure of 2 x 10 ^ Torr, 

_ 7 

although it was capable of attaining vacuums in the middle 10 Torr 
range. A general photograph of our entire equipment setup is shown 
in figure 3. 

The major components of the vacuum system are one Eimac HV-1 
diffusion Pump, one fore pump, two liquid air cold traps, an 80 cm. 
octoil manometer, a high vacuum Alpert valve, a glass manifold, and two 
ionization gauges. 

The system generally performed in a satisfactory manner; it was a 
relatively tight system and quite easy to keep clean; however, all 
operations were greatly slowed down by the Eimac diffusion pump. The 
operating base pressure of 2 x 10 ^ Torr was a compromise between a 
desired high degree of purity, and the very low pumping rate of the 
diffusion pump, which was operated with the baffle assembly installed. 
This mode of operation reduced the pumping speed by a factor of two, 
but it prevented pump oil from backing up into the rest of the system. 
Although this type diffusion pump is theoretically capable of attaining 

_7 

a vacuum of 5 x 10 Torr, the pumping speed at this pressure is less 
than 20/' of the pumping speed at 2 x 10~^ Torr. 

Pump performance was improved by substituting Dow Corning 704 

silicone fluid for Eimac Type A diffusion pump fluid. A pressure of 

—6 

2 x 10” Torr could be attained in about five to seven days using the 
following procedure. 
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Figure 2 VACUUM SYSTEM SCHEMATIC 
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EjjgttEe General Equipment' Layout 



The discharge tube, and portions of the manifold were wrapped with 
electro-thermal Model H-T 362 heating tapes, initially set at 50 volts. 
When equilibrium between the pumping rate and the outgassing rate was 
reached, the voltage on the heating tapes was increased. This was re- 
peated ’until a maximum volta. e of 125 volts was reached. The system 
was allowed to bake in this configuration for about 24 hours. The 
electrodes were cleaned by heating to a red-orange color with an in- 
duction heater. All filaments were outgassed by heating them with 
the conventional filament power supply unit, (The power supplies are 
discussed in more detail in the subsequent section of this paper , ) 

Both of these operations were repeated until equilibrium between the 
pumping and outgassing rates \tfas established. The heating tapes were 
operated simultaneously with the electrode heating operation to min- 
imize deposites of outgassed materials being formed on glass portions 
of the system. Final clean up was attained by operating an Argon dis- 
charge, then evacuating the system. This final step was repeated until 
the Argon discharge attained its characteristic violet glow. 

The system's pressure is measured using an ionization gauge VG-1A 
triodi . Argon pressure in the discharge tube was measured using the 
Octoil S manometer. (One cm. oil equals 0,672 mm, Hg,). 

The U-shaped discharge tube was made of Vycor glass. The electrodes 
were fabricated from Tungsten filament attached to Tantalum supports; 
circular anodes of Tantalum were similarly mounted. The electrode as- 
semblies were mounted on VG-1A triode bases,, which were fused to the dis- 
charge tube. This mode of mounting electrodes is very inflexible; cathode 
replacement is a very inefficient and time consuming process. Discussion 
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of our experiments with other schemes for mounting the electrodes 
will be found in section 3 of this paper. 

2.2 i-lectrical System. 

The general circuit layout is shown in figure 4. Three power 
supplies were used. A Kepco, Model 770~B, 600 volt 3 amp. unit provided 
the high voltage potential for the discharge tube. Fine discharge 
current control was obtained using a resistor bank. 

The filaments were heated with a Kepco, Model KM 236=! 5A 50 volt, 

15 amp. unit. No external control devices were used with this power 
supply. 

A Kepco, Model 500R-B 600 volt, 0.3 amp. unit was used to provide 
regulated power for the Pupps or auxiliary discharge maintained at the 
anode end of the discharge tube. A sliding wire resistor was used in 
this circuit to facilitate current control of the Pupps discharge. 

2.3 Electronic System. 

Detection of moving striations can be accomplished with a rotating 
mirror or by electronic methods. We used both systems; the rotating 
mirror gave us gross information relating to the presence of the moving 
striations and their stability. 

Detailed measurements of various striation parameters such as 
frequency and wavelength were made using the electronic circuit shown 
schematically in figure 5« The major components were a Tektronix type 
54 5 A oscilloscope, two RCA type 7102 multiplier phototubes, and a battery 
pack. 

The RCA-7102 tube is well-suited for use in the near~infrared region 
that was of interest to us. Maximum response of the RCA-7120 is the region 
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Figure 4* Electrical System 




Figure 5» Electronic System 
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‘igure 5* i hototube Assembly 



0 + 0 

8000A - 1000A , our field of interest fell in this v;avelength intena L <> 

A battery pack was used to produce a 1200 volt DC source for the photo 
tubes. Figure 6 shows the method used to mount the tubes. 

2.4 Spectrographic Equipment. 

In order to be able to measure a Doppler shift in an Argon gas 
discharge, a Fabry-Perot interferometer was mounted on the optical b ;nch 
of a Bausch Lomb, Littrow type spectrograph. The general arrangement 
is shown in figure 7. 

The spectrographs dispersion readily separated the interferometric 

patterns when entrance slit width was set at 30 microns. V/ider entrance 

o o 

slit settings caused overlapping of the Argon lines at 8006A and 80-5A. 

All exposures were made on type I-N Kodak spectrographic plates because 

o 

of their sensitivity in the 800QA region of the spectrum. A nominal ex- 
posure time with this arrangement, an Argon discharge, and a discharge 
current of 300 ma was three minutes. 

The Fabry=P 3 rot interferometer was constructed for maximum resolution 
in the 7500A to 10,0Q0A region of the spectrum. The details of its 
construction are contained in Appendix I. The method of calibration and 
the results are contained in section 5» 
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Feasibility of using a r e movable end plate tube for zas discharges <> 

3.1 Problem. 

The technique by which gas discharge tubes had been fabricated was 
reviewed prior to mounting the "U" tube on the vacuum system for the 
Doppler experiment. Up to this time VGIA ion gauge bases had served as 
mountings for the electrodes, these bases being fused to pyrex tubing 
of desired length and diameter tc form the closed discharge tube. This 
technique had distinct disadvantages. For example, the tube had to be 
completely disassembled when the cathode burned out, a common occuranee 
caused by ion bombardment at high currents. Also different electrode 
arrangements and materials could not be tried out because of the limited 
amount of tine available. In order to overcome these limitations an 
attempt was made to establish an argon gas discharge in a tube with re= 
movable end plates. 

3.2 Procedures and Results. 

In an a ttempt to increase cathode lifetime a "L" dispenser cathode 
with approximate dimensions of 3b x 36 x 12 mm had been obtained. Pre- 
vious experiments indicated that an anode configuration as illustrated 
in figure 8 might induce stability to the amplitude of moving striations 
as observed on a photomultiplier. As a result of these two considerations, 
a 40 mm diameter tube was fabricated to an overall length of 92.3 err. using 
commercial Pyrex 0-ring joints and pyrex tubing. In order to silver solder 
"Hermetic Seal" lead=in pins and facilitate maching, end plates were made 
of brass. The completely assembled tube is illustrated in figure 9. 

The assembly was first tested on a Veeco vacuum leak detector. The 
tube was evacuated to the 10 ^ Torn range in about ten minutes. Spray- 
ing the exterior of the end plates and seals with helium indicated neg= 



f 5 
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Figure 8. Anode Configuration , 
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Figure 9. Removable Endplate Assembly 
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ligible leakage,, The tube was then fused to the vacuum system* evacuated, 
and baked out with heating tapes,, A base pressure of 5 x 10 Torr was 
obtained. 

An attempt was made to activate the dispenser cathode; however, this 
attempt failed primarily due to inexperience with activation procedures 
(i.e, decomposition of strontium carbonates was too fast due to ex- 
cessive filament temperature), A helical wound tungsten filament with 
a cylindrical shield was substituted for the dispenser cathode. 

With the system re-evacuated and baked out to a base pressure of 
2 x 1CT 5 Torr, the tube was filled with argon to a pressure of 5 mm Hg, 
When a discharge could not be obtained using a 600 volt power supply, 
a 1200 volt, 0,75 ma, DC power supply was placed across the tube, A 
discharge was obtained with this power supply; however, after about a 
minute of operating a break-down occurred between the electrode and 
the end plate at the cathode which exceeded the current capacity of 
the power supply. The end plates were removed and the anode plate 
was coated by brush with glyptol insulating varnish and the cathode 
plate with a dissolved lucite dope. The tube was reactivated; however, 
the lowest obtainable base pressure was 3 x 10 Torr, The tube was 
filled with argon, but on this occasion the discharge became unstable 
and extinguished itself with a change of color that indicated possible 
contamination. It is relieved that leakage was occur ing around the "0" 
rings due to the rough nature of the dielectric coatings. At this point, 
because of ti e Imitations, the removable end plate investigation was 
discontinued and VGIA bases were used to assemble electrodes with thor- 
iated tungsten filaments. These bases were then fused onto the n U" tube. 



J? 



Using a center intense image the spacer may be ground and adjusted with 
spring tension to parallelism within about a hundredth of a wavelength 
of the reference light. Thus any further changes of the fringe pattern 
as the eye is moved across the field would be due to lack of optical 
flatness of the plates. 

When the spacer was ground as exactly as possible the aluminized 
plates were removed and the dielectric coated plates were substituted 
in their place. The interferometer was set up an the optical bench 
of a Bausch and Lomb Littrow type spectrograph. 

The interferometer was adjusted using a Neon Geissler tube and 
ground glass diffuser as a light source. Sufficient red light in the 
visible spectrum was reflected by the dielectric coating so as to be 
produce a weak interference pattern. The eye when focused at infinity 
could distinguish this pattern against the background intensity suff- 
iciently well to allow adjustment of the interferometer springs for 
parallism. 

Next a white light source was placed at the plate holder of the 
spectrograph in order to project the optical axis from the entrance 
slit of the spectrograph. The focusing lens was positioned by aut©= 
collimaticn using this source and the image of the spectrograph slit 
projected onto itself. 

The next problem incurred was that of placing to interferometer 
on its table such that the vertical diameter of interference pattern 
was coincident with the slit of the spectograph. Reflectance of the 
interferometer plates in the visible region was insufficient to allow 
projection of the interferometric pattern onto the entrance slit. 



Although the removable end plate tube did become operational, it 
is felt that this method of assembly is practical. 

3.3 Conclusions 

It is recommended that the following steps be taken, (ij A 25 mm 
tube be assembled of 90 cm length in order to reduce the current in the 
discharge. (2) End plates be made of 5 " or 3/4 phenolic insulating 
board rather than brass. (3) Lead-in seals be set in with vacuum epoxy 
using 4 parts 828 epoxy resin and 1 part Shell Epon curing agent "U" . 
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4. Investigation into Two Methods for Obtaining Radial Distribution 
of Emitters from Spectroscopic Measurements. 

4.1 Problem. 

It is possible to investigate certain parameters of optically-thin 
plasmas using spectroscopic techniques. These methods leave the plasma 
virtually undisturbed; however, the data are frequently very difficult 
to interpret. Barr^^ and Bockasten^^ have reported on methods to 
convert observed radiation intensity profiles from an optically thin 
cylindrical source into an expression for the radial distribution of 
emitters. Basically, both methods are readily adaptable to computer 
programs for the Control Data Corporation Model 1604 Computer at 
USNPS. 

From the nature of the mathematical transformations that were used 
in the development of both systems, there appears some question regarding 
the sensitivity of the observed intensities to small changes in the 
radial distribution. There was also a need for curves of expected in- 
tensities from assumed radial distribution. These two questions could 
be answered simultaneously, with the results yielding some insight into 
the relative accuracy of each method. 

4.2 Procedure. 

A highly colimated detector mounted as shown in figure 10, accepts 
th energy from a very narrow tube of light through the source . 
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Figure 10 

Plasma Disk showing coordinate system 
(The 2 axis is directed into the paper) 

The radial distribution of emitters, E(r), is assumed to be circularly 
Sjf.iuneti ic with respect to the £ axisj the radiation is assumed to be 

isotropic, and not absorbed by the plasma. Under these conditions, 

Jo 

XM Axaz =<S^E(r) AxAu (4.1) 

'-f=o 

KJ 

Expression (4.1) reduces to the integral equation (4.2) when infinitely 
small volume elements are considered. 
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!(*)=£ I E(r)dy 



(4.2) 



2 2 2 

Using the fact that y = r - x , equation (4.2) becomes: 

-R 



X('%)= 2. 



ECr) rdr 



(4.3) 



Using Abel's transformation this becomes: 

, R 

E ( r ) = 



H. 



dx> 



(4.4) 



ft" / r 

The observations that yield l(x) in practice do not give a smooth 
curve or an easily plotted function, but only a series of numerical 
values. There are various numerical methods to convert these empirical 
values into a smooth curve, l(x). Since equation (4.4) uses the d erivative 
of I(x). This system possesses certain inherent insensitivities which 
must be investigated. 

Equation (4.3) can be used to evaluate expected observed intensities 
from assumed radial distributions. On can evaluate l(x) exactly from a 

. i* 2 

well behaved function such as E(r) = 1 - ( ZA . ) j however, such simple 

J 

functions are too few in number to be of much help; therefore numerical 
methods must be used. 

A relatively simple expression for l(x) can be obtained by fitting 
E(r) to a second degree polynominal in the same manner as Bar7 

2,0 

mmm* 

Tk 



79 / 



2,0 

I(*k) "A 21 a J* E(fi) 



(4.5) 



where A equals .05 and the coefficients A are determined according 

jK 

to (4.6). 



a jk. ~ */$ ( 2 k + ) ) 



'A 



if 
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A comparison between the values of l(x) for E(r) ■ i-K) 2 ob- 
tained from equations (4.3) and (4.5) is shown in figure 11. Figures 
12, 13 > and 14 show expected l(x) curves determined by equation (4.5) 
from 15 assumed radial distributions. In all graphs, both E(r) and 
I(x) are normalized so that their maximum values are equal to unity. 

All transformations are based on 20 equal increments ( ), with E(r) 

reducing to zero at the outer boundary of the plasma, both l(x) and 
E(r) having zero slopes at the origin. The numerical values used to 
plot figures 12, 13, and 14 are tabulated in table 2 of Appendix II. 

We did not have a satisfactory discharge tube for making intensity 
observations to compare the Barr and Bockasten methods. Therefore the 15 
sets of calculated values for l(x) were used as inputs to determine the 
corresponding radial distribution using both methods. The computed values 
are tabulated in table 3 of Appendix II. 

4.3 Results and Conclusions. 

The errors inherent to the numerical methods employed in equation 4.5 
do not appear excessive. The maximum error shown figure 11 is 1.5$, tut 
the average value at each calculated point is 0.1761$. Therefore, we 
assumed that the numerical methods used in converting the assumed radial 
densities into profiles are at least as accurate as the experimental 
measurements. 

The 1 5 radial distributions were not selected randomly, but with 
three objectives. First, some were used because they closely approximated 
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Figure t2 
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physically measured distributions. The zero order Bessel function (curve 
1) and the Gaussian distribution (curve 2) are examples of this type. 
Second, families of similar curves were selected in order to determine the 
sensitivity of intensity curves to changes in radial distribution. Curves 
7-11 and 12-15 are examples of this method of selection. Finally, some 
"limiting value" distributions were selected. They are not physically 
possible, but their discontinuities do represent the boundary conditions 
within which all distributions must lie. Curves 6 and 13 are examples 
of this type. 

It is apparent from an inspection of figures 12, 13, and 14 that 
the observed intensity profile does not change nearly as much as the 
radial distribution function. Note that in the family of curves 7-11, 

E(0) varies from 0 to 0.8, but l(0) varies only from 0.8 to 1.0. One 
may conclude that the density of emitters at the center of the cylinder 
is not a maximum value if the intensity profile is not maximized at 
x = 0. The converse is not true, 1(0) may be the maximum value of l(x), 
and a rarefaction may occur at r = o. (curve 11). 

A general, qualitative picture of the radial density may be obtained 
by curve fitting an experimental plot of intensities to those curves 
obtained in this paper from the assumed radial distributions; however, 
quantitative values of electron or atom density can not be obtained by 
this method. 

A comparison between the methods of Barr and Bockasten using these 
assumed curves is somewhat artificial, since the curves are not all phy- 
sical] y possible, and the discontinuities strain both methods beyond 
their designed capabilities. An inspection of table 3 of Appendix II 
shows that neither method is very reliable near the outer boundary of 
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the region containing the emitters. Generally, the errors produced by 
the Barr method are larger than those produced by Bockasten's method by 
a factor of at least two. The errors in both cases change sign frequently, 
indicating an oscillation or crossing of the exact curve by the generated 
curve. 

Barr's met .od incorporates a least squares correction to improve the 
quality of input data. This causes an anticipation of large changes in 
radial distributions, and hence a rounding off of the generated intensity 
curves. Consequently Barr's errors are larger on a point -by-point, com- 
parison, but the overall curves are just as valid as Bockasten's curves. 

The inversion of the l(x), profiles can give, at least, only re- 
lative values of the radial variation of the emitters, E(rj. The inter- 
pretation of the radial density of emitters in terms of the radial 
variation of electron density, for instance, requires certain assumptions 
about the uniformity of temperature, density of atoms, etc. The mag- 
nitude of the electron density must be obtained by some other method, 
such as measurements of absolute intensity. 
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5. Doppler Shift Experiment 

5.1 Problem. 

The problem for which this experiment was designed is discussed 
in Section 1 .4. To review briefly, a Doppler shift should exist if 
the phenomena of moving striations in gas discharges is due to the 
collective motion of the light sources. Therefore, the objective was 
to make a measurement of sufficient precision to determine if such a 
Doppler shift occurs. 

The experimental lay-out consisted of a "U" shaped gas discharge 
tube with windows at the base which provide a means to view the pos- 
itive column. On a line with this section of the discharge was the 
optical bench of a spectrograph upon which was mounted the necessary 
focusing lens and a Fabry-Perot interferometer. A shift in the inter- 
ferometric pattern when the striations were moving in one direction 
versus moving in the other direction would indicate the Doppler shift, 

5.2 Instrumentation Parameters. 

The initial parameters for the experiment were established in the 
following manner. From the work of A. W. COOPER it appeared feasible 
to obtain moving striations in a gas discharge in Argon with a velocity 
of 90 meter/sec. By comparing striations moving toward to those moving 
away from the interferometer a difference of velocity of 180 m/sec. could 
be obtained. As previously mentioned, the work of Robertson and Hakeem 
established the dependence of striations on the metastable state. Accord- 
ingly as part of the experimental procedure the wave lengths in Argon 
upon which the measurement would be attempted would include decays to the 
metastable states. In A -gon such wavelengths are in the near red region 
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of the spectrum. The expected doppler shift would be: 



A A = (f)A = ai,s* o.oosa 

For a Fabry-Perot interferometer the intensity, I, at any point 
in the field is the following function of I q , the uniform monochromatic 
source intensity, 



1= L 



i 

I + q 2 - s/n 2 d/^ 



where 



Cl 2 - 






0 -R )*• 
d = 2TT (~^A cos 6^) 



(5.0 



R is the reflectance of the plates, t is the thickness of the spacer, 
is the index of refraction of the gas in the space between the 
plates^ A is the wave length of light. The half width of an inter- 
ference fringe is given by quanity Q S I n i ) . For a high 

resolution instrument this value may be obtained as a small perturbation 
of an order of interference. 

d=2TT(n + fi) (5.2) 

Where n is an integer^ £ is a perturbation 

Sin 4^ = sinrrn cosne 1- cosrrn s/nTfs 

— TT d for £ sufficiently small 



Therefore ^4 

a its - ± 1 (5.3) 

But , by taking the differential of the argument 



£ = 




^ ~ S 4 n cos 9 

; ' a 2 



dA 



(5. a ) 



£ = - n 



dA 

A 
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Substituting 5.3 into 5.4. 




A £ 

£ t jj na 



Assuming that the index of refraction is unity, that plates have 99% 

o 

reflectance and a wavelength of interest is 8115A the minimum spacer 

o 

thickness for measuring a 0.005A shift is: 



A & __ ( 6//SA) 

£ 7T Q 4A (£ ty )(£oq)('o,oo5A) 



In order that spacer thickness would not be the limit of resolution the 
latter was made approximately a order of magnitude larger than the cal- 
culated minimum. 

The instrumental broadening caused by lack of flatness of the plates 
can be calculated to a first approximation. The dielectric plates were 
specified flat to within 1/100 of a wave length of a Hg green line. The 
resulting change in order of interference Sm of the plates due to lack 
of flatness St is: m A — 

A “ £ <5t 



S 





(5.5) 



where Ao is wavelength of Hg green line. The expected change of the 
order of interference due to Doppler shift Am is: 



A A 

~ = X~ 

Thusjthe ratio of the changes in order is: 



8m _ a / 


' Ac \ 


Arm rn AA 1 


too ) 


— / Ao \ 


c 


l 100 ) 


i v 


5 Goo A 


/ 3 X 1 0 s 'Vsec A 


10 0 


\/ cm Xl80 7 5 ec7 



(5.6) 



(5.7) 
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For a one centimeter spacer, the change of order of interference due 
to irregularities in the flatness of the plates amounted to approximately 
eight tenths the order of interference change due to Doppler shift. 

Thus, the expected instrumental broadening was O.CXJ+A. Obviously this 
condition placed a lower limit on practical spacer size. 

Dielectric coated plates with the specified flatness were obtained 
from Perkin Elmer Corp. and a 8.679 i 0.001 mm spacer cvas manufactured 
of Invar. The interferometer was assembled as described in Appendix I, 

5.3 Interferometer Calibration. 

Zeeman splitting of A-gon was used to determine the minimum resolution 

of the Fabry-Perot interferometer. An Argon Geissler tube was located 

normal to the axis of a six inch circular pole piece^g-inch gap electro- 

magnet. The magnet provided cross fields of up to three kilo gauss. An 

image of the source was focused onto the optical axis of the interferometer 

and spectrograph. Exposures were taken at various values of magnetic 

field strength. Line traces were then made on the interferometric 

pattern us:ing a knorr-Alber micro photometer. With a 25mm slit height 

only two interference fringes were imaged on the photographic plate. 

Expected Zeeman spitting patterns for lines in the I'egion of interest are 

shown in figure 15. These patterns were calculated for a nominal field 

/ 1 2/ 

strength of 10,000 gauss. ' Several line patterns were analyzed and 

o 

the best resolution was obtained on argon 801 5 A where the splitting could 

be observed at 300 gauss. This corresponds to a wavelength shift of 
o 

0.01 5 A, a factor of three larger than the desired resolution. 

Certain problems became apparent while making the calibration plates, 

o 

First, line widths were approximately 0.03A. It was not determined at 
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Figure 15 
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this point whether the lines were instrumental broadened or source 
broadened. Second, it was necessary to stop the microphotometer down to 
a small aperture, approxi- atley 0e5mm x 0.05mm in order separate adjacent 
lines. To compensate for the small aperture the sensitivity of the instru- 
ment was raised to the maximum value. As a result the emulsion "noise" level 
showed significantly on the line trace. The results of the calibration runs 
indicate that at least a factor of three must be gained in the resolution of 
the interferometer. 

5.A Results and conclusions. 

An attempt was made to complete the experiment; however, other factors 
degraded the results. Several casualties were sustained in the discharge 
tube while trying to conduct the experiment. On only one occassion was the 
experiment completed satisfactorily. However , on that occassi on the striaticn 
velocity was 26 m/sec, a factor of four below the desired 90 n/sec. 4 velo- 
city '±^ high as 52 m/sec, was obtained during one of the unsuccessful runs. 

Ho Doppler shift could be detected in the 26 m/sec. run. 

An overall factor of eight must be gained in the system performance 
before satisfactory results nay be anticipated. In order to obtain this im- 
provement, the answer must be determined as to v:hether the broad line width 
of trie interferometric patterns is caused by source broadening or instrumental 
broadening. The interferometer should be checked against a narrow line 
width source such as a Schuler tube to make this determination. If the broad 
width is caused by the source, the interferometer resolution may be improved 
by increasing spacer size.. 

liven with improved interferometer resolution, the measurement of small 
■o pier shifts on broad lines imposes problems. These migh* be 
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solved by using a narrow line width absorbsion cell in the optical path 
or by using t power points of off center images in the Fabry-Perot field. 
This latter method consists of adjusting the interferometer so that the 
line of interest just appears as a first interference ring. The "tail" 
of the line intersects the axis of the interferometer at the £ power 
point as shown in figure 16. The height of the center portion of the 
pattern would then be a very sensitive measurement of the broadness 
of the line and would reflect a Doppler shift should it exist. 

An increase of signal to noise ratio on the microphotometer traces 
is another improvement which must be made in either case. This might 
be obtained oy using an emulsion with fine grain, such a type IV plate . 




Figure 16 

quarter Power Point Resolution of Fabry-Perot Interferometer 
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6. Recommendations 



It is recommended that the following areas be investigated by 

future research groups at USNPGS* These recommendations are eased 

upon problems we encountered and subjects of interest that became 

apparent during the course of our experimenting* 

The removable endplate concept should be reviewed, and investigation 

made into locating an endplate material that would perform satisfactorily 

“7 

as an insulator and maintain a vacuum of 10 Torr* 

Experiments should be made to compare the spectroscopic method of 
obtaining emitter distribution with probe and/or microwave techniques* 
Spectroscopic data should be reduced by noth the Barr and Bcckasten 
methods in order to gain further insight regarding the relative merits 
of each method 

Further attempts should be made to determine if an orderly motion 
of atoms occurs when moving striations are observed in a glow discharge. 
The discharge tube used for this experiment should have a smaller inside 
diameter than the tube we used, because a smaller tube will produce moving 
striations with a greater velocity, thus producing a greater Doppler 
effect* 



37 



7 , Acknowledgements 

The authors wish to acknowledge the assistance of Professor Ro L 0 
Kelly for his counsel and advice, and Professor A. Vf c Cooper for his aid 
in maintenance of the vacuum system® We are also indebted to Professor 
Jo Pc Duff in for his assistance in preparing the computer programs 0 

We also gratefully acknowledge the technical assistance of Messers 
r.o Smith, R« C. Moeller, and J # CaLder. 



38 



Bibliography 

1. Sanborn C. Brown, Speech, Sigma Xi - KESA National lecture, I960 



2. J. J. Thompson and G. P, Thompson, Conduction of Electricity 
Thro ug h Gases . (Cambridge University Press, 1933 J 3rd. Ed., Vol II. 

3. A. W. Cooper, Moving S triations in the Inert Gases . (Ph.D, Thesis, 
Queens University of Belfast, 1961) 

4. Sanborn C. Brown, Basic Data of Plasma Physics . (The Technology 
Press and Wiley and Sons, Inc., New York, 1959), Chap. 14. 

5. F. H. II. Kinley and R. J. O'Malia, A Spectroscopic Inves tigat ion 
o f An A rgo n Glow Discharge . (M. o. Thesis, U. S. Naval Postgraduate 
School, 19*62)”." 

6. Gorder: Francis, The Glow Discharge at Low Pressure . (Handbuch der 
Physik, Springer - Verlag, 1956), Vol XXII. 

7. H. 5. Robertson and K. A. Hakeem, Proceedings of the Fifth Inter- 
natio na l Conference on Ionization P he nomena in Gases . (North Holland 
Publishing Company, Amsterdam, 1962) p. 550. 

8. U. Asccli - Bartoli, A, De Angelis, and S. Martellucci, Proceedi ng s 
of the F ifth International Conference on Ionization Phenomena in 
Gases . (North - Holland Publishing Company, Amsterdam, 1962 ) P. 376. 

9. W. L. Barr, J. Opt. Soc . An. 52, 885 (1962). 

10. K. Bockasten, J. Opt. Soc, Am. 51, 943 ( 1 96 1 ) . 

11. H. Margenau and G. M. Murphy, The Mathematics of Physics and Chemistry ,, 
(D. Van Ho strand Company, Inc., Princeton, New Jersey, 19431 P> 507. \ 

\ 

12. J. C. Van Den Bosch, The Zeeman Effect . (Handbuch der Physik, Springer 
Verlag, 1956), Vol. XXVIII, p.“296. 

13. °. Tolansky, High Resolution Spectroscopy . (Pitman Publishing Corp,, 
New York, 1945;, p. ”138.*” 

f 

i 

? 



39 



* 




Appendix I 

Manufacturing and Adjusting a Fabry-Perot Interferometer 
for Use in Near Red Spectrum 

In order to observe certain near red Argon line shifts a Fabry- 

o 

Perot interferometer was constructed to be used in the 7500 - 1 0 , 000 A 
region. Two 2-inch diameter fused silica plates with dielectric coat- 
ings were obtained from the Perkin-Elmer Corp. for the etalon. The 
transmittance of these plates are shown in figure 17. Various problems 
were incurred in manufacture and adjustment of the interferometer be- 
cause of the very small reflectance in the visible spectrum. The most 
expeditious and sensitive solutions of these problems are reported here 
for future reference. 

The housing for the etalon was fabricated in the physics department 
shop and is illustrated in figure 18. Two 2-inch diameter glass optical 
flats were coated with aluminum thin films by vacuum vaporization tech- 
niques to an overall reflectance of approximately 50 %. A spacer was 
made by inserting three Invar pegs into a 2-inch diameter retaining ring 
at peripheral points corresponding to the pressure points of the inter- 
ferometer housing. The pegs were ground on a surface plate until they 
were equal within 0.001" as measured with a micrometer. The inter- 
ferometer was assembled with the aluminized plates and spacer. A low 
pressure mercury lamp was used as a light source to determine the need 
for further spacer grinding in order that the plates would be optically 

parallel. This determination was made by conventional methods and is 

/ 1 3/ 

outlined here . 

From the fringe pattern, the high point of the spacer can be 



ik 0 



Figure I? 

Reflectance of Plates Used in Fabry-Perot Interferometer 




U1 




c 






b2 



rv'r^ 



determined and the spacer subsequently ground to reduce the wedge angle. 
Figure 19 (a) shows a moderate degree of non -parallelism part (b) shows 
a reduce degree of non-parallelism, and part (c) shows nearly parallel 
plate pattern. 




Fart (a) Fart (b) Part (c) 

Figure 19 

Course Adjustment of Fabry-Perot Interferometer 
The final grinding is determined by moving the eye across the field of 
the interferometer. If the center most fringe opens when the eye is 
moved toward one peripheral point and closes when moved to the dia- 
metrically opposite point, the point toward which the fringe opens is 
the high point (see figure 20). 





Figure 20 

Fine Adjustment of Fabry-Perot Interferometer 
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However, a sensitive method of accomplishing this alignment was that of 
using the white light source at the plate holder and observing multiple 
reflections of the slit back on to the slit holder. Figure 21 illustrates 
this adjustment. 




Interferometer Interferometer 

axis non parallel Aligned axis non parallel 

Figure 21 

Centering of Fabry-Perot Pattern on Spectrograph 
The interferometer was adjusted on the table so that the center ref- 
lection was coincident with the entrance slit of the spectrograph. 

A hard piece of "tacky wax" was placed around each foot of the inter- 
ferometer prior to aligning. Once aligned, the wax was heated with 
a heat gun and then allowed to cool. This proved to be an excellent 
way of securing the interferometer against accidental jarring and at 
the same time did not disturb the interferometer alignment in the pro- 
cess. 
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APPENDIX II 



TABLE 1 

COMPARISON BETWEEN EXACT SOLUTION AND NUMERICAL METHOD SOLUTION 





I(x) 


I(x) 


ERROR 




FROM EQ (5.3) 


FROM EQ (5.5) 




1 


1.000000 


1 .000000 


0,0000 


2 


0.99248 


0.99628 


+0.3828 


3 


0.98496 


0.98506 


+0.0101 


4 


0,96241 


0.96647 


+0.4218 


5 


0.93984 


0.94063 


+0.0840 


6 


0.90977 


0.90775 


-0.2220 


7 


0.86466 


0.86811 


+0.3990 


8 


0.83458 


0.82202 


-1.5049 


9 


0.76692 


0.76989 


+0.3872 


10 


0.71428 


0.71221 


-0.2898 


11 


0.65413 


0.64954 


-0.8270 


12 


0.57969 


0.58254 


+0.4916 


13 


0.5H27 


0.51202 


+0.1467 


14 


0.43985 


0.43887 


-0.0182 


15 


0.36842 


0.36422 


-1.1400 


16 


0.28947 


0.28939 


-0.0276 


17 


0.21804 


0.21600 


-0.9356 


18 


0.14586 


0. 14618 


+0.2194 


19 


0.08271 


0.08282 


+0.1329 


20 


.03082 


0.03044 


-1.2337 


21 


0.00000 


0.00000 


0.0000 
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TABLE 2 



CALCULATED INTENSITIES FROM ASSUMED RADIAL DISTRIBUTION 





CURVE 1 


CURVE 2 


CURVE 3 


CURVE 4 


CURVE 5 


.i.K 


E(r) 


I(x) 


E(r) 


I(x) 


E(r) 


I (xl. 


E_(rJ _ 


I(x) 


E(r) 


I(x) 


1 


1 .0000 


0.9933 


1.0000 


0.6287 


1.0000 


1.7146 


1.0000 


1.3333 


1 .0000 


1.0021 


2 


0.9604 


0.9828 


0.9802 


0.6163 


1.0000 


1.7117 


0.9975 


1 .3283 


0.9500 


0.9909 


3 


0.9120 


0.9575 


0.9233 


0.5807 


1.0000 


1.7029 


0.9900 


1.3134 


0,9000 


0.9662 


4 


0.3463 


0.9218 


0.8353 


0.5257 


1.0000 


1.6881 


0.9775 


1.2886 


0.8500 


0.9315 


5 


0.8075 


0.8816 


0.7262 


0.4573 


1.0000 


1.6672 


0.9600 


1.2541 


0.8000 


0.8890 


6 


0.7652 


0.8336 


0.6067 


0.3823 


1.0000 


1 .6399 


0.9375 


1.2103 


0.7500 


0.8401 


7 


0.7196 


0.7788 


0.4868 


0.3068 


1.0000 


1.6059 


0.9100 


1.1574 


0.7000 


0.7860 


8 


0.6711 


0.7183 


0.3700 


0.2365 


1 .0000 


1.5647 


0.8775 


1.0960 


0.6500 


0.7278 


9 


0.6201 


0.6532 


0.2780 


0.1767 


1.0000 


1.5159 


0.S400 


1.0265 


0.6000 


0,6664 


10 


0.5669 


0.5847 


0.1980 


0.1263 


1.0000 


1 .4585 


0.7975 


0.9496 


0.5500 


0.6027 


11 


0.5112 


0.5140 


0.1354 


0.0869 


1.0000 


1.3915 


0.7500 


0.8660 


0.5000 


0.5373 


12 


0.4554 


0.4423 


0.0889 


0.0580 


1.0000 


1.3134 


0.6975 


0.7767 


0.4500 


0.4711 


13 


0.3980 


0.3707 


0.0612 


0.0374 


1.0000 


1.2219 


0.6400 


0.6827 


0.4000 


0.4049 


14 


0.3400 


0.3005 


0.0341 


0.0217 


1.0000 


1.1136 


0.5775 


0.5851 


0.3500 


0.3395 


15 


0,2818 


0.2330 


0.0198 


0.0126 


1 .0000 


0.9818 


0.5100 


0.4856 


0.3000 


0.2756 


16 


0.2239 


0.1698 


0.0110 


0.0070 


0.9900 


0.8138 


0.4375 


0.3858 


0.2500 


0.2143 


17 


0. 1 666 


0.1123 


0.0060 


0.0037 


0.9500 


0.5805 


0.3600 


0.2880 


0.2000 


0.1566 


18 


0.1104 


0.0626 


0.0030 


0.0018 


0.6000 


0.2702 


0.2775 


0.1949 


0.1500 


0.1038 


19 


0.0555 


0.0234 


0.0015 


0.0009 


0.1000 


0.0409 


0.1900 


0.1104 


0. 1000 


0.0577 


20 


0.0025 


0.0010 


0.0008 


0.0003 


0.0010 


0.0004 


0.0975 


0.0406 


0.0500 


0.0208 


21 


0.0000 


0.0000 


0.0000 


0.0000 


0,0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 
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APPENDIX II 



TABLE 2 (CONTINUED) 





CURVE 6 


CURVE 7 


CURVE 8 


CURVE 9 


CURVE 10 


J,K 


E(r) 


I(x) 


E(r) 


I(x) 


E(r) 


I(x) 


E(r) 


I(x) 


E(r) 


I (x) 


1 


0.0000 


0.9635 


0.0000 


0.8964 


0.2000 


0.9458 


0.4000 


0.9887 


0.6000 


1.0602 


2 


0.0000 


0.9672 


0.0200 


0.9125 


0.2100 


0.9565 


0.4100 


0.9936 


0.6050 


1 .0623 


3 


0.0000 


0.9789 


0.1700 


0.9632 


0.3000 


0.9917 


0.4200 


1 .0096 


0.6350 


1.0698 


4 


0.0000 


0.9996 


0.4350 


1 .0346 


0.5000 


1.0461 


0.5000 


1.0461 


0.7000 


1 ,0813 


5 


0.0000 


1.0317 


0.8000 


0.0917 


0.8000 


1.0917 


0.8000 


1.0917 


0.8000 


0.0917 


6 


0.0000 


1 .0799 


0.9500 


1.0867 


0.9500 


1 .0867 


0.9500 


1.0867 


0.9500 


1.0867 


7 


0.0000 


1.1557 


1.0000 


1.0420 


1 .0000 


1.0420 


1 .0000 


1.0420 


1.0000 


1.0420 


8 


0.0000 


1.3065 


0.9950 


0.9698 


0.9950 


0.9698 


0.9950 


0.9698 


0,9950 


0.9698 


9 


1 .0000 


1.5159 


0.9800 


0.8755 


0.9800 


0.8755 


0.9800 


0.8755 


0.9800 


0,8755 


10 


1 .0000 


1.4585 


0.9150 


0.7533 


0.9150 


0.7533 


0.9150 


0.7533 


0.9150 


0.7533 


11 


1 .0000 


1.3915 


0.8000 


0.6087 


0.8000 


0.6087 


0.8000 


0.6087 


0.8000 


0.6087 


12 


1.0000 


1.3134 


0.6000 


0.4576 


0.6000 


0.4576 


0.6000 


0.4576 


O.60OO 


0,4576 


13 


1 .0000 


1.2219 


0.4450 


0.3357 


0.4450 


0.3357 


0.4450 


0.3357 


0,4450 


0.3357 


14 


1.0000 


1.1136 


0.3200! 


0.2383 


0.3200 


0.2383 


0.3200 1 


0.2383 


0.3200 


0,2383 


15 


1 .0000: 


0.9818 


0.2300 


0.1627 


0.2300 


0.1627 


>0.2300 


0.1627 


0.2300 


0.1627 


16 


0.9900 


0,8138 


0.1500 


0.1032 


0.1500 


0.1032 


0.1500 


0.1032 


0.1500 


0.1032 


17 


0.9500* 


0.5805 


0. 1000 


0.0608 


0. 1000 


0.0608 


0.1000 


0 . 0608 


0.1000 


0.0608 


18 


0.6000 


0.2702 


0.0500 


0.0287 


0.0500' 


0.0287| 


; 0.0500 


0.0287 


0.0500 


0.0287 ; 


19 


0.1000 


0.0409 , 


■0,0200 


0.0115 


0.0200 


0.0115 


0.0200 


0.0115 


0.0200 


0,0115 


20 


0.0010 


0.0004 


0.0100 


0.0042 


0.0100 


0.0042 


0.0100 


0.0042 


0.0100 


0,0042 


21 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


O.OOOo 


0.0000 


O.OuOO 


0.0000 
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TABLE 2 (CONTINUED) 





CURVE 11 


CURVE 12 


CURVE 13 


CURVE 14 


CURVE 15 


J,K 


E(r) 


I(x) 


E(r) 


I(x) 


E(r) 


I(x) 


E(r) 


Kx) 


E(r) 


I(x) 


1 


0.8000 


1.1257 


1 .0000 


1.0447 


0. 1000 


1.9504 


1 .0000 


1.3159 


1 .0000 


1.6174 


2 


0.8050 


1.1242 


1.0000 


1 .0395 


1 .0000 


1.9479 


1 .0000 


1.3119 


1.0000 


1 .6143 


3 


0.8100 


1.1192 


1.0000 


1.0236 


1.0000 


1.9401 


1 .0000 


1 .2999 


1.0000 


1.6048 


4 


0.8250 


1.1114 


1 .0000 


0.9964 


1.0000 


1 .9272 


1.0000 


1 .2795 


1 .0600 


1.5890 


5 


0.8500 


1.1017 


1 .0000 


0.9563 


1 .0000 


1.9090 


1 ,0000 


1 .2504 


1.0000 


1.5664 


6 


0.9500 


1.0867 


1 .0000 


0.8997 


1 .0000 


1 .8852 


1 .OoOO 


1.2116 


1 .0000 


1.5370 ! 


7 


1.0000 


1.0420 


0.9700 


0.8202 


1.0000 


1 .8558 


1 .0000 


1.1620 


1 ,0000 


1.5001 


8 


0.9950 


0.9698 


0.8900 


0.7189 


1 .0000 


1 .8204 


1 .0000 


1 ,0992 


1 .0000 


1 .4553 


9 


0.9800 


0.8755 


0.7850 


0.6048 


1 .0000 


1.7788 


1.0000 


1.0181 


1,0000 


1,4016 


10 


0.9150 


0.7533 


0.6500 


0.4840 


1 .0000 


1.7303 


0.9600 


0,9H5 


1,0000 


1.3378 


1 1 


0.8000 


0.6087 


0.5000 


0.3682 


1 .0000 


1.6745 


0.8800 


0.7841 


1 .0000 


1 .2622 


12 


0.6000 


0.4576 


O.36OO 


0.2696 


1 .0000 


1 .6105 


0.7700 


0.6434 


1 ,0000 


1,1718 


13 


0.A450 


0.3357 


0.2500 


0.1942 


1 .0000 


1.5374 


0.6200 


0.4978 


1 .0000 


1,0602 i 

j 


14 


0.3200 


0.2383 


0.1800 


0.1402 


1 .0000 


1.4538 


0.A550 


0.3667 


0,9700! 


0.9166 j 


15 


0.2300 


1 

0.1627 


0. 1300 


0.0989 


1 .0000 


1 .3577 


0.3300 


0.2651 


0.8750 1 


0.7396 


16 


0.1500 


0.1032 


0.0900 


O.O607 


1.0000 


1.2462 


0.2400 


0.1872 


0 

CO 

i> 

0 

0 


0.5448 


17 


0. 1000 


0.0608 


0.0600 


0.0429 


1 .0000 


1.1146 


0.1700 


0. 1264 


0.5100 j 


0.3583 


18 


0.0500 


0.0287 


0.0400 


0.0257 


1.0000 


0.9547 


0,1200 


0,0793 


0.3300 


0.2103 


19 


0.0200 


0.0115 


0.0230 


0.0127 


1 .0000 


0.7478 


0.0720 


0.0412 


0.1860 


0.1028 


20 


0.0100 


0.OJ42 


0.0100 


0.C042 


1.0000 


0.4163 


0.0350 


0.0146 


0.0800 


0.0333 


21 


0.0000 


0 .0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0,0000 
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TABLE 3 

COMPARISON BETWEEN RADIAL DISTRIBUTION DETERMINED BY BARR AND BOCKASTEN METHODS 







CURVE 1 






j 


E(r) 


E(r)* 


ERROR* 


E(r)** 


ERROR** 


1 


1 .0000 


0.9595 


-4.0470 


1.0015 


0. 1470 


2 


0.9604 


0.9483 


-1.2568 


0.9632 


0.2926 


3 


0.9120 


0.9153 


0.3607 


0.9124 


0.0482 


4 


0.8463 


0.8625 


1 .9201 


0.8494 


0.3639 


5 


0.8075 


0.8144 


0.8533 


0.8081 


0.0706 


6 


0.7652 


0.7683 


0.4038 


0.7657 


0,0601 


7 


0.7196 


0.7230 


0.4683 


0.7200 


0.0597 


8 


0.671 1 


0.6737 


0.3859 


0.6715 


0.0596 


9 


0.6201 


0.6232 


0.5089 


0.6205 


0.0645 


10 


0.5669 


0.5693 


0.4322 


0.5839 


2.9970 


1 1 


0.51 12 


0.5143 


0.6084 


0.5117 


0.0958 


12 


0.4554 


0.4578 


0.5292 


0.4557 


0.0703 


13 


0.3980 


0.4007 


0.6859 


0.3984 


0,1005 


14 


0.3400 


0.3429 


0.8588 


0.3404 


0. 1206 


15 


0.2818 


0.2846 


1 .0043 


0.2822 


0.1526 


16 


0.2239 


0.2269 


1 .3354 


0.2244 


0.2144 


17 


0.1666 


0.1698 


1 .9088 


0.1672 


0.3361 


18 


0.1104 


0.1142 


3* 4783 


0.1113 


0.7880 


19 


0.0555 


0.0643 


15.9279 


0.0552 


-0,5225 


20 


0.0025 


0.0276 


- 


0.0095 


- 


21 


0.0000 


0.0071 




0.0000 


0.0000 
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CURVE 2 






j 


E(r) 


E(r)* 


ERROR* 


E(r)** 


ERROR** 


1 


1.0000 


0.9905 


-0.9540 


1.0010 


0.0980 


2 


0.9802 


0.9733 


-0.6968 


0,9813 


0,1102 


3 


0.9233 


0.9225 


-0.0856 


0.9243 


0,1083 


4 


0.8353 


0.8391 


0.4621 


0.8366 


0.1563 


5 


0.7262 


0.7355 


1,2723 


0.7278 


0,2162 


6 


0.6067 


0.6201 


2.2105 


0,6087 


0.3412 


7 


0.4868 


0.5030 


3,3215 


0.4883 


0,3122 


8 


0.3700 


0.3925 


6.0728 


0.3731 


0,8297 


9 


0.2780 


0.2954 


6.2768 


0,2797 


0.6043 


10 


0.1980 


0.2144 


8.2862 


0.2016 j 


r 

1,7875 


11 


0.1354 


0.1505 


11,1915 


0.1368 


1,0490 


12 


0.0890 


0.1020 


14.5971 


0.0908 


2,1014 


13 


0.0612 


0.0667 


9,0746 


0.0612 


0,1308 


14 


0.0341 


0.0418 


22,7633 


0,0352 


3.1974 


15 


0.0198 


0.0249 


25.6061 


0.0203 


2.3232 


16 


0.01 10 


0.0140 


26.7452 


0.0113 


2,9918 


17 


0.0060 


0.0077 


28.7854 


0.0062 


2.9950 


13 


0.0030 


0 . 004 1 


37.3333 


0.0032 


6,6667 


19 


0.0015 


0.0021 


38,6667 


0.0014 


-3.3333 


20 


0.0008 


0.0009 


20.6667 


0,0007 


-606667 


21 


0.0000 


0.0004 


<=, 


0.0000 


0,0000 



* BARE 



** BOCiOiSTEN 



50 



APPENDIX II 



TABLE 3 (CONTINUED) 







CURVE 3 






j 


E(r) 


E(r)* 


ERROR* 


E(r)** 


ERROR** 


1 


1.0000 


1.0018 


0.1800 


1 .0000 


0,0000 


2 


1 .0000 


1 .0009 


0.0870 


1,0000 


0,0000 


3 


1 ,0000 


0.9996 


-0.0370 


1.0000 


O.UOOO 


4 


1.0000 


0.9976 


-0.2380 


1 ,0000 


0,0000 


C 

-/ 


1 .0000 


1 .0000 


0.0000 


1 ,0000 


0,0000 


6 


1.0000 


1.0003 


0.0300 


1.0000 


0.0000 


7 


1,0000 


1.0012 


0.1200 


1,0000 


0.0000 


8 


1.0000 


0.9995 


-0.0440 


1.0000 


0,0000 


9 


1.0000 


1 .0016 


0. 16o0 


1.0000 


0,0000 


10 


1 .0000 


1 .0005 


0,0460 


1,0550 


5.4950 


1 1 


1 .0000 


1.0001 


0.0160 


0.9999 


-0.0080 


12 


1.0000 


0.9996 


-0.0440 


0.9998 


-0.0160 


13 


1.0000 


1 .0000 


| 0.0000 


0.9996 


-0,0360 


14 


1.0000 


0.9987 


| -0,1310 


0,9994 


-0.0060 


15 


1.0000 


0.9879 


-1.2090 


0,9975 


-0.2470 


16 


0.9900 


0.9395 


-5.0990 


0,9935 


0.3545 


17 


0.9500 


0.7981 


-15.9916 


0.9323 


-1,8652 


18 


0,6000 


0.5569 


-7.1833 


0,5833 


-2,7833 


19 


0.1000 


0,2924 


192,4100 


0.1327 


32,7300 


20 


0.0010 


0.0959 


- 


0,0147 


- 


21 


0.0000 


0,0120 


— 


0.0000 


0,0000 
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CURVE 4 






j 


E(r) 


E(r)* 


ERROR* 


E(r)«* 


EKkGR** 


1 


1 .0000 


1.0008 


0.0840 


1 

0,9999 


-0.0090 


2 


0.9975 


0.9978 


0.0271 


0.9975 


0,0000 


3 


0.9900 


0.9894 


-0.0626 


0.9900 


0,0000 


4 


0.9975 


0.9758 


-0.1964 


0.9975 


0,0000 


5 


0.9600 


0.9596 


-0.0365 


0,9600 


0,0000 


6 


0.9375 


0.9372 


-0.0267 


0.9375 


0,0000 


7 


0.9100 


0,9105 


0,0516 


0.9100 


0,0000 


8 


0,8775 


0.8769 


-0.0627 


0,8775 


0,0000 


9 


0.8400 


0.8407 


0.0881 


0.8400 


0.0000 


10 


0.7975 


0.7973 


-0,0188 


0,8282 


3 8501 ' 


11 


0.7500 


0,7498 


-0.0267 


0.7500 


0,0000 


12 


0.6975 


0.6970 


-0.0674 


0,6975 


0,0000 


13 


0,6400 


0,6399 


-0.0156 


0,6400 


0.0000 


14 


0.5775 


0.5775 


0,0000 


0.5775 


0,0000 


15 


0,5100 


0,5095 


-0.0863 


0.5100 


0.0000 


16 


0.4375 


0.4372 


-0.0594 


0.4375 


0,0000 


17 


0,3600 


0,3598 


-0,0500 


0.3600 


0.0000 


18 


0.2775 


0,2775 


0.0000 


0,2775 


0.0252 


19 


0.1900 


0,1912 


0.6421 


0.1905 


0 . 2684 


20 


0,0975 


0.1090 


11.8256 


0,0936 


-4,0308 


21 


0.0000 


0.0467 




0,0000 1 


0,0000 
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CURVE 5 








r E(r) 


E(r)* 


ERROR * 


I E(r)** 


ERROR** 


1 


1.0000 


0.9490 


-5.1040 


1.0012 


0.1190 


2 


0,9500 


0.9384 


-1 .2158 


0,9521 


0,2210 


3 


0,9000 


0.9075 


0.8367 


0.9010 


0,1177 


4 


0.8500 


0.8571 


0.8412 


0,8504 


0.0529 


5 


0.8000 


0.8070 


0.8762 


0.8001 


0.0137 


6 


0.7500 


0.7560 


0.7960 


0.7498 


-0,0307 


7 


0.7000 


0.7056 


0.8057 


0.6994 


-0,0900 


8 


0.6500 


0.6543 


0.6662 


0,6488 


=0.1815 


9 


0,6000 


i 0.6046 


0.7683 


0.5992 


=0.1333 


10 


0.5500 


0.5536 


0.6582 


0.5497 


=0,0491 


1! 


0.5000 


0.5033 


0.6580 


0.5250 


5,0000 


12 


0.4500 


0.4528 


0.6333 


0.4531 


0,6867 


13 


0.4000 


0.0428 


0.7075 


0.4003 


0 . 0800 


14 


0.3500 


0.3527 


0.7600 


0.3503 


C.0857 


15 


0.3000 


0.3022 


0.7467 


0.3003 


0.0900 


16 


0,2500 


0.2522 


0.8760 


0.2502 


0.1000 


17 


0.2000 


0.2021 


1 .0500 


0.2002 


0. 1200 


18 


0.1500 


0,1521 


1 ,3867 


0.1502 


0.1667 


19 


0.1000 


0.1026 


2.50OO 


0.1004 


0,4500 


20 


0.0500 


0.0574 


14.9800 


0.0483 

d 


=3,3000 


21 


0,0000 


0.0243 


=. 


0.0000 i 


0.0000 



* BARR 

** BOCKASTEN 



53 



APPENDIX II 



TABLE 3 (CONTINUED) 







CURVE 6 




j 


j 


L(r) 


L(r)* j 


r 

| ERROR* 


E(r)** 


ERROR'"* 

» - ... ■■ - J 


1 


OoOOOO 


0.0012 j 


1 1 


-0,0028 


1 


2 


0.0000 


0.0006 




0,0023 


- | 




0.0000 


-0.000 1 


- 


0.C005 


. 1 


4 


0.0000 


-0,0019 


! 


0.U006 


- 


s 


0.0000 


0,0004 




C.C014 


- 


6 


0.0000 


0.C095 


1 


0.0063 


- 


7 

i 


0.0000 


C.C858 


_ 


0,00273 


; 


8 


0.0000 


0.3173 


- 


O.C895 




Q 

/ 


1.0000 


0.6215 


-37.8480 


C.9174 ! 


-8 0 2640 * 


1C 


1 .0000 


0.88JC 


-1 1.9C40 


1 .0550 

I 


5.4950 j 


1 1 1 


1 .0000 


J .0001 


0.0160 


o 

o 

o 

o 

o 


I 0.0000 j 


! 12 


1 .0000 


i 

0.9996 


) 

=0.0440 


1 

0.9998 | 


j 1 

j -0.0160 j 


i 13 


1 .0000 


1 .0000 1 


0,0000 


0,9996 | 


1 -0,0360 j 


! 14 


1 .0000 


0,9987 


-0.1310 1 
l i 


j 0,9994 


J {j 

I -0.C580 j 


| 15 


1 .0000 


0.9879 


i 

-S.2090 


1 0,9975 


j =0.2470 J 


\ 

16 


0.9900 


0.9395 


1 1 

j =5.0990 


I * 

j 0,9935 


0,3545 1 

& 


17 


0.9500 ; 


0.7981 


-15.9916 


0.9323 


-1,8652 i 

J 


18 


0,6000 


0.5569 


=7.1833 

! 


0.5833 


-2. 7833 | 


19 


0,1000 ! 


0.2924 


j 192,4100 


0.1327 


| 32.730oJ 


i 20 


0.C01C 


0.0959 




O.C147 


! * ! 


i 

1 21 


0.0000 


0.C120 


i 


0.0000 


j i 

j, 0.0000J 
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CURVE 7 


j 


E ( r ) 


E ( r )* 


ERROR * 


E ( r )** 


ERROR * 


1 


0.0000 


- 0.0155 


- 


- 0,0023 


- 


2 


0.0200 


0.0339 


69.7000 


0.0178 


- 1 1 .0500 


3 


0.1700 


0.1819 


7.0000 


0.1648 


- 3.0529 


A 


0,4350 


0.4196 


“ 3.5402 


0.4366 


0,3586 


5 


0,8000 


0.6746 


= 15.6750 


0.781 1 


- 2 . 36 o 2 


6 


0.9500 


0.8636 


- 9.0958 


0.9418 


- 0.8631 


7 


1 .0000 


0.9601 


- 3.9900 


0.9951 


- 0,4890 


8 


0.9950 


0.9790 


- 1.6030 


0,9947 


- 0.0271 




0.9800 


0.9535 


- 2,7010 


0,9769 


= 0,3133 


10 


0.9150 


0.8838 


- 3.4098 


0,9291 


1.5387 


1 1 


0.8000 


0.7700 


= 3.7413 


0,7941 


- 0,7400 


12 


0.6000 


0,6249 


4.1500 


0.6047 


0,7833 


13 


0 . 4450 


0,4760 


6,9753 


0.4479 


0.6427 


14 


0.3200 


0.3469 


8,4062 


0.3236 


1 , 118 ? 


15 


0,2300 


0,2461 


6,9913 


0,2308 


0,3478 

! 


16 


0.1500 


0, 1683 


12,2133 


0,1530 


1 , 9867 


17 


0.1000 


0.1081 


8.1500 


0.1000 


0.0000 


18 


0.0500 


0.0632 


26.4600 


0,0514 


2,9200 


19 


0,0200 


0,326 


63.0000 


0.0218 


8.8500 


20 


0,0100 


0.0138 


38,4000 


0.0097 


- 3,4000 


21 


0,0000 


0,0048 


^ 


0.0000 


0,0000 
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— . — — — — — — — . — — — — — — 

CURVE 8 


J 


E(r) 


E(r)* 


ERROR* 


E ( r )** 


ERROR** 


1 


0.2000 


0.1668 


-16.1100 


0.1992 


-0.3850 


2 


0.2100 


0.2052 


-2.2714 


0.2084 


-0.7476 


3 


0.3000 


0.3173 


5.7667 


0.2974 


-0,8633 


4 


0.5000 


0.4985 


-0,2860 


0.5020 


0.3920 


5 


0.8000 


0.7070 


-11.6275 


0.7864 


-1.7050 


6 


0.9500 


0.8711 


-8.3011 


0.9417 


-0.8737 


7 


1.000 


0.9601 


-3.9900 


0,9951 


-0,4900 


8 


C.9950 


0.9790 


-1.6030 


0.9947 


i -0,0271 


9 


0.9800 


0.9535 


-2.7010 


0.9769 


-0,3133 


10 


0.9150 


0.8838 


-3.4098 


0.9291 


1.5388 


11 


0.8000 


0.7700 


-3.7413 


0.7941 


-0.7400 


12 


0.6000 


0.6249 


4.1500 


0.6047 


0.7833 


13 


0.4450 


0.4760 


6.9753 


0,4479 


0.6427 


14 


0.3200 


0.3469 


8.4000 


0.3236 


1.1187 


15 


0.2300 


0.2461 


6.9913 


0.2308 


0.3478 


16 


0.1500 


0.1683 


12.2133 


0.1530 


1.9867 


17 


0.1000 


0.108? 


8.1500 


0,1000 


0.0000 


18 


0.0500 


0.0632 


26.4600 


0.0514 


2.9200 


19 


0.0200 


0.0326 


63.OOOO 


0.0218 


8,8500 


20 


0.0100 


0.0138 


38.4000 


0.0097 


-3.4000 


21 


0.0000 


0.0048 


- 


0.0000 


0.0000 
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CURVE 9 


0 


E(r) 


E(r)* 


ERROR* 


E(r)** 


ERROR** 


1 


0.4000 


0.3484 


=12.8900 


1 

0,4009 


0.2275 


2 


0.4100 


0.3709 


-9.5268 


I 0.4103 


0.0829 


3 


0.4200 


0.4382 


4.3357 


| 0,4164 


-0,8548 


4 


0.5000 


0.5519 


10.3840 


0.51 14 


2.2740 


5 


0 . 8000 


0.7203 


-9.9588 


0,7863 


=1,7150 


6 


0.9500 


0.871 1 


-8.301 1 


0.9418 


1 =0,8631 


7 


1.0000 


0.9601 


-3.9900 i 


! 0.9951 1 


| -0.4890 j 


8 


0.9950 


0.9790 


-1.6030 


0.9947 1 


=0.0271 


9 


0.9800 


0.9535 


=2,7010 


0,9769 


-0.3 133 


10 


0.9150 


0.8838 


=3.4098 


0.9290 


1 .5388 


1 1 


0.8000 


0.7700 


-3.7413 


0.7941 


-0.7400 


12 


0.6000 


0.6249 


4J500 


0.6047 


0.7833 


13 


0.4450 


0.4760 


6,9753 


0.4479 


0.6427 


14 


0.3200 


0.3469 


8.4000 


0,3236 


1.1187 


15 


0.2300 


0.2461 


6,9913 


0,2308 


0.3478 


16 


0.1500 


0.1683 


12.2133 


0.1530 


1.986? 


17 


0.1000 


0. 1081 


8.1510 


0.1000 


0.0000 


18 


0.0500 


0.C632 


26.4600 


0.0515 


2.9200 


19 


0.0200 


0.0326 


63.0000 


0.0218 


8.8000 


20 


0.0100 


0.0138 


38.4000 


0.0097 


-3.4000 


21 


0.0000 


0.0048 




0.0000 1 


0.0O00 
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CURVE 10 






3 


| E(r) 


E(r)* 


_ 

ERROR* 


E(r)** 


ERROR** 


1 


0.6000 


0.5888 


=1.8700 


0.5996 


“6,7317 


2 


0.6060 


0.6016 


-0.5570 


0.6044 


=0.0926 | 


3 


0.6350 


0.6401 


0.8016 


0.6346 


=0.0567 


4 


0.7000 


0.7048 


0.6843 


0.6980 


=0,2886 


5 


0.8000 


0.7994 


-0.0763 


0.8021 


0.2675 


6 


0.9500 


0.8944 1 


j =5.8516 


0,9416 


=0.8821 


7 


1 ,0000 


0.9601 


=3 0 9900 


0,9949 


=0,5080 


8 


0.9950 


0.9790 


=1 ,6030 


0,9945 


=0,0492 


9 


0 9800 


0.9535 


-2.7010 


0 9767 


=0.3377 


10 


0.9150 


0,3838 


-3 ,4098 


0,9288 


1.5082 


1 1 


0,8000 


0.7700 


-3 7413 


0,7937 


-0.7812 i 


12 


0.6000 


0.6249 


4J500 


0.6043 


0,7150 


13 


0,4450 


0.4760 


6,9753 


0,4473 


0,5213 


14 


0,3200 


0.3469 


8.4000 


0,3228 


0.3750 


15 


0.2300 


0.2461 


6.9913 


0,2302 


0,0913 


16 


0,1500 


0.1683 


12.2133 


0, 1438 


-4.I4OO 


17 


0.1000 


0.1081 


8.1500 


0. 11 23 


12.3400 


13 


0.0500 


0.0632 


26.4600 


0.0528 


5,6800 


19 


0,0200 


0.0326 


63,0000 


0.0218 


8,8500 


20 


0,0100 


0.0138 


38,4000 


0 0097 


=3 .4000 


21 


0.0000 






- . ,&flQ0Q 
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CURVE 11 


J 


E(r) 


E(r)* 


ERROR* 


E(r ) 


ERROR** 


1 


0.8000 


0.7985 


-0.1813 


0.8000 


0.0000 


2 


0.8050 


0.8018 


-0.4025 


0,8048 


-Q.C285 


3 


0.8100 


0.8014 


0.1778 


0.8105 


0.0654 


4 


0.8250 


0.8294 


0.5394 


0.8230 


-0,2364 


5 


0.8500 


0.8700 


2.4529 


0.8543 


0.5178 


6 


0.9200 


0.9225 


-2,8989 


0,9458 


=0.4358 


7 


1.0000 


0.9660 


-3=3980 


0.9951 


-0,4890 


3 


0.9950 


0.9790 


-1 .6030 


0,9947 


-0,0271 


9 


0.9800 


0.9535 


-2,7010 


0.9769 


1 =0,3133 


10 


0.9150 


0.8838 


-3.4098 


I 0.9291 


1.5388 


11 


0.8000 


0.77o0 


=3.7413 


0.7941 


-0.7400 


12 


0.6000 


0.6249 


4J500 


0.6047 


0.7833 j 


13 


0.4450 


0.4760 


6.9753 | 


0,4479 


0.6427 


14 


0.3200 


0.3469 


8.4000 


0,3236 


1,1187 


1 15 


0.2300 


0.2461 


6,9913 


0.2308 


0.3478 


16 


0.1500 


0.1683 


12.2133 


0,1530 


1.9867 


17 


0.1000 


0.108! 


8.1500 


0.1000 


0,0000 


18 


0.0500 


0 .0632 


26.4600 


0,0514 


2,9200 


19 


0.0200 


0.0326 


63.0000 


0.0218 


8.8500 


20 


0.0100 


0.0138 


38.40000 


0,0097 


-3.4000 


21 


0.0000 


0.0048 1 


_ - 


0,0000 


0,0000 
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CURVE 12 


j 


E(r) 


E(r)* 


ERROR* 


I E(r)** 


ERROR** 


1 


1.0000 


1.0008 


0.0800 


0.9998 


-0.0230 


2 


1 .0000 


1 .0005 


0.0500 


0,9998 


-0.0220 


3 


1 .0000 


-0.0190 


0.9998 


-0.0240 


-0.0470 


4 


1 .0000 


0.9989 


=0.1130 


0,9995 


-0.0470 


5 


1.0000 


0.9972 


-0.2790 


1.0000 


0,0000 


6 


1.0000 


0.9837 


-1.6340 


! 0.9986 


-0.1360 


7 


0.9700 


0.9473 


-2.3422 


0.9672 


-0.2897 


8 


0.8900 


0.8777 


=1.3798 


0.8897 


-0.0371 


9 


0.7850 


0.7768 


l 

-1.0433 


1 

0.7844 


-0.0802 


10 


0.6500 


O.65IO 


0.1523 


0.O590 


1 .3892 


11 


0.5000 


0.5160 


3J940 


0.5018 


0,3600 


12 


0.3600 


0.3871 


7.5222 


0.3628 


0,7694 


13 


0.2500 


0.2795 


11.7960 


0.2536 


1 ,4280 


14 


0.1800 


0,1984 


10.2222 


0.1819 


1.061 1 


15 


0.1300 


0. 1403 


7.9615 


0.1310 


0,7538 


16 


0.0900 ' 


0,0985 


9.5000 


0.0908 


0,9333 j 


17 


O.ObOO j 


0.0673 | 


12.1667 


O.O0O9 


1.5333 


18 


0.0400 


0.0438 1 


9,5750 


0.0403 


0,7250 


19 


0.0230 


0,0263 


14.1739 


0.0234 


1.7391 


20 


0.0 100 | 


0.0134 


34.2000 


0.0101 


0,8000 


21 


1 

0.0000 i 


0.0052 


- _ 


0,0000 


0,0000] 
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CURVE 13 


j 


E(r) 


E(r)* 


ERROR* 


E(r )** 


ERROR* 


1 


1 .0000 


1 .0019 


0.1900 


0.9994 


-0.0560 


2 


1.0000 


1 .0009 


0.0900 


0.9994 


-O.O62O1 


3 


1.0000 


0.9995 


-0,0500 


0.9993 


=0.0650 ! 


4 


1 .0000 


0.9975 


-0.2500 


0.9993 


-0,0660j 


5 


1 .0000 


1.0000 


0,0000 


0,9993 


=0,0650 


6 


1 ,0000 


1.0003 


0.0300 


0.9993 


-0.0700 ' 


7 


1 .0000 


p 

1,0012 


0.1270 


0.9993 


-0,0740 


8 


1 .0000 


j 0.9994 


-0,0600 


0.9992 


I -0.0770 


9 


1 .0000 


1 ,0025 


0.2460 


0.9992 1 


-0,0830 ; 


10 


1.0000 


1.0007 


0.0740 


1,0683 


6,8270 j 


11 ; 


1.0000 


1 . 0003 


0.0310 


0,9990 


=0, 0980 | 


12 


1 .0000 


0.9996 


-0.0360 


0.9989 


-0. f 130 


13 


1 ,0000 


1 ,0007 


0.0670 


0.9987 


=0,1310 


14 


1 .0000 


1 . 0007 


0.0720 


0.9984 


=0,1550 


15 


1 .0000 


0.9996 


-0.0370 


0,9980 


=0,1970 


16 


1 .0000 


0,9998 | 


-0,0200 


0.9973 


-0,2650 


17 


1.0000 


0,9995 


-0,0520 


0.9959 


=0,4120 


18 


1 .0000 

I 


0,9883 


-1 , 1670 


0,9898 


-1.0210 : 


19 


1.0000 


0,9001 


=9,9900 


1.0183 


1 ,8270 j 


20 


1 ,0000 


0.6633 


=33,6690 


0,8269 


=17,3080 


21 


0,0000 


0,3680 


- 


0 . 00 jO 


0.0000 1 



* BARR 



** BOCKASTEM 



6 J 



APPENDIX II 



TABLE 3 (CONTINUED) 



CURVE 14 



j 


E(r) 


| E(r)* 


ERROR* 


E(r)** 


| ERROR*-* 


1 


1 .0000 


1,001 1 


0.1110 


0.9995 


-0,0050 ! 


2 


1.0000 


1.0007 


0.0670 


0.9998 


-0.0140 


3 


1.0000 


0.9998 


-0.0230 


0,9999 


0.0080 


4 


1 .0000 


0.9985 


-0. 1530 


0.9999 


-0,0060 


5 


1.0000 


0.9999 


-0.0050 


0,9998 


-0,0150 


6 


1 .0000 


0.9999 


-0.0090 


0.9998 


-0.0190 


7 


1 .0000 


1 ,0003 


0,0270 


0.9995 


-0.0530 


8 


1.0000 

1 


! 0.9953 


-0.4660 


1.0003 


o.02eo 


9 


i.ooo) ! 


| 0.9795 j 


-2.0500 


0,9795 


0,0000 


10 


0.9600 1 


I 0.9373 | 


-2.3677 


0.9801 


2,0896 


1 1 


0,8800 


0.8631 


-1.9239 


0.8792 


-0.0954 


12 


0.7700 


0,7543 


-2.0403 


0,7685 


-0.1948 


13 | 


0.6200 


0.6218 


0.2935 


0.6195 


-0.0855 


14 j 


0.4550 


0.4835 


6,2571 


0.4587 


0.8088 j 


15 


0.3300 


0.3585 


8.6333 


0.3332 


1 0.9818 


16 1 


0.2400 


0.2593 


8.0292 


0.2417 


0.7167 


17 j 


0.1700 


0. 1846 


8.5941 


0, 1720 


1.1647 


18 | 


0.1200 


0.1266 


5 o 4833 


0, 1202 


0.1833 


19 


0,0720 


0.0797 


10.6806 


0,0732 


1.6944 


20 


0.0350 


0.0425 


21 .3714 


0,0341 


-2.5142 


JLL 


0.0000 


0.0172 




0,0000 


0.0000 1 
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1 CURVE 


15 


j 


E(r) 


E(r)* 


ERROR* 


E(r)** 


ERROR** 


1 


1 .0000 


1,0016 


0,1610 


0.9999 


r 

-0.0130 j 


! 2 


1 .0000 


) .0008 j 


0.0810 ' 


0.9999 j 


! I 

-0,0030 


3 


1 . 0000 


0,9997 


-0.0330 


1.0000 


0,0000 


1 4 


1 ,0000 


0.9978 


-0,2230 


0,9999 


-0,0030 


5 

j 


1 .0000 


0.9999 ; 


’ -0,0050 


0,9999 


-0.0030 1 


6 


1.0000 


1.0002 


0,0190 


0,9999 


-0.0030 


7 


1 .0000 


1.0010 


0,0130 


0.9999 


-0.0030 


a 


1 .0000 


0,9996 


-0.0420 


0.9999 


-0,0030 


9 


1 .0000 


1 .0014 


0.1420 


! 0.9999 


-0,0030 


10 


1.0000 j 


1 ,0003 


0.0270 

i 


1.0475 

1 


4,7480 


ii 


1 .0000 


0.9996 


-0.0350 

t 


| 0,9995 
| 


-0,0500 


J 12 
t 


1.0000 


j 0.9955 


-0,4440 


0.9998 


-0.0160 : 


; 13 


1,0000 1 


| 0.9791 


-2 . O 890 


0=9984 


! -o , 1 640 1 

1 


j 14 


0,9700 


0.9322 


-3,9000 


0,9660 ; 


! -0,4082; 


15 


0,8750 


0.8387 


=4, I486 


0.8721 


-0.3360 


16 


0,7180 


0.6985 


-2.7159 


0,7147 


-0.4512 


17 


0.5100 


0.5298 


3.8804 


0,5130 


0.5843 


18 


0.3300 


0.3614 


9.5121 


0.3333 


0,9970 


i 19 


0.1860 


0.2168 


16.5430 


0,1896 


1.9624 


! 20 
t 


0.0800 


0. 1090 

j 


36,2875 


0=0810 


1.2375 


! 

LJ-L. 


I 0.0000 


O.C419 




1 0,0000 


0.0000 
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